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INVERSE  SCATTERING  APPLICATION 
IN  DETERMINING  TERRAIN 
FEATURE  PARAMETERS 


INTRODUCTION 

Purpose.  The  purpose  of  this  research  note  is  to  present  the  results  of  a  pre¬ 
liminary  study  of  inverse  scattering  and  its  application  in  determining  terrain  feature 
parameters. 


Background.  The  direct  scattering  problem  consists  of  solving  for  a  reflected 
or  scattered  wave,  given  a  wave  incident  on  a  particular  scatterer.  In  the  solution  to 
this  problem,  all  the  necessary  physical,  electrical,  and  statistical  properties  of  the 
scatterer  are  assumed  to  be  known,  which  then  can  be  used  to  construct  a  model  to 
calculate  the  scattered  fields.  Sometimes,  however,  various  properties  of  a  scatterer 
are  not  known.  An  example  of  this  is  the  simulation  of  a  layer  of  vegetation  with  a 
continuous  random  medium.  To  compute  the  scattered  electric  field  from  the  random 
medium,  one  must  know  the  two-point  correlation  function  of  the  dielectric  fluctua¬ 
tions.  This  correlation  function  has  never  been  determined  for  vegetation  of  any  type, 
and  so  the  difficulty  of  calculating  the  scattered  fields  remains. 


The  inverse  scattering  problem  is  defined  as  using  various  properties  of  the 
scattered  fields  to  compute  certain  physical  or  electrical  properties  of  a  scatterer. 
Several  examples  of  inverse  scattering  problems  will  be  considered  briefly  in  this 
section.  Later,  in  the  analysis  section,  some  of  these  examples  will  be  developed  in 
more  detail. 


One  example  of  an  inverse  scattering  problem  is  calculating  the  physical  profile 
characteristics  of  a  discrete  object  of  revolution  using  the  far  field  measurements  of 
the  radar  cross  section.  In  this  problem,  a  discrete  object  of  revolution  is  said  to  exist, 
and  one  makes  radar  scattering  measurements  as  a  function  of  incidence  angle  on  the 
object.  The  radar  scattering  measurements  are  used  to  compute  the  radar  cross  section, 
which  in  turn  is  used  to  calculate  the  profile  of  the  object  via  a  particular  set  of  in¬ 
version  equations.  Two  basic  techniques  can  be  used  to  solve  this  problem,  geometric 
optics  and  physical  optics.  The  objects  of  interest  are  usually  assumed  to  have  in¬ 
finite  conductivity,  which  is  a  good  assumption  if  one  is  trying  to  determine  the 
profile  characteristics  of  an  aircraft,  missile,  or  satellite. 
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Another  example  of  an  inverse  problem  is  calculating  the  dielectric  constant 
of  a  slightly  rough  surface  using  measurements  of  the  radar  backscatter  coefficient  for 
both  horizontal  and  vertical  polarizations,  along  with  the  first-order  pertubation 
scattering  models  for  a  slightly  rough  surface.  Because  the  dielectric  constant  of  a  soil 
can  be  related  to  moisture  content,  a  solution  to  this  inverse  problem  presents  the 
possibility  of  determining  quantitatively  the  soil  moisture  content  in  vast  areas  of 
desert  type  terrain. 


Another  example  of  an  inverse  problem  is  calculating  the  dielectric  constant  and 
conductivity  of  materials  located  beneath  a  forest  canopy.  A  bistatic  radar  setup  is 
used  to  measure  amplitude  and  phase  of  the  reflected  field.  The  measurements  are 
then  used  in  the  appropriate  inversion  equations  to  calculate  the  conductivity  and 
dielectric  constant  of  the  medium  below  the  canopy.  These  electrical  properties  can 
then  be  used  to  give  a  good  indication  of  the  nature  of  the  material  that  is  present. 
The  derivation  of  these  inversion  equations,  along  with  a  complete  analysis  of  this 
problem,  will  be  presented  later. 


Another  type  of  inversion  problem  involves  solving  a  Fredholm  linear  integral 
equation  of  the  first  kind.  Such  integral  equations  arise  in  determining  dielectric 
constant  profiles  as  a  function  of  soil  depth  or  in  computing  the  temperature  profile 
of  the  atmosphere.  The  solution  of  these  integral  equations  involves  reducing  the 
integral  equation  to  a  system  of  linear  algebraic  equations  that  must  be  solved  by 
matrix  inversion  techniques.  A  great  deal  of  research  has  been  conducted  on  the 
solution  of  the  Fredholm  integral  equation  and  its  difficulties,  such  as  stability  and 
uniqueness. 


In  the  analysis  section,  more  detailed  discussions  will  be  presented  of  the  in¬ 
version  problems  mentioned  above.  Computational  results  will  be  presented  for  the 
slightly  rough  surface  problem  and  the  forest  canopy  problem.  In  each  of  these  two 
cases,  the  geometry  of  the  problem  to  be  solved  will  be  presented,  and  the  appro¬ 
priate  inversion  equations  will  be  derived.  A  discussion  section  will  follow,  in  which 
side-looking  airborne  radar  imagery  will  be  considered  for  extracting  terrain  features. 
This  is  another  example  of  an  inversion  problem  as  applied  to  terrain  feature  deter¬ 
minations. 
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An  infinitely  conducting  object  is  formed  by  revolving  the  profile  r(z)  about 
the  z  axis.  A  radar  at  point  A  illuminates  the  object  and  measures  the  backscattered 
power  in  such  a  way  that  the  radar  cross  section  can  be  determined  for  the  look 
direction  77.  The  objective  of  this  problem  is  to  determine  the  object  profile  r(z) 
from  the  remote  sensing  measurements  of  radar  cross  section  as  a  function  of  tj. 
Bojarski  has  derived  a  solution  to  this  problem  using  a  geometric  optics  approach 
that  consists  of  the  following  two  parametric  equations: 1 


r2(i?)  = 


(17')  sin  tj'  cos  rj'  dr)' 


z(v) 


1  f  q(7?>) 

*  J  r(rj') 


sin2  17'  dr)' 


(1) 


(2) 


where  o(rj)  is  the  radar  cross  section  as  a  function  of  look  direction.  The  advantage 
of  using  the  above  equations  is  that  no  phase  information  is  required  to  obtain  the 
profile  r(z).  The  disadvantage  is  that  the  radar  cross  section  must  be  known  for  all 
values  of  17  from  0  to  n.  In  addition,  there  are  restrictions  on  the  curvature  of 
r(z)  in  order  for  the  above  equations  to  hold. 


Bojarski  has  also  derived  expressions  to  determine  r(z)  for  an  object  using  a 
physical  optics  approach.2  The  starting  point  for  this  derivation  is  the  equation  for 
the  scattered  magnetic  field  given  below . 


H 

—  S 


(B  x  U) 


x  V0  ds 


(3) 


'sorbert  N.  Bojarski.  Signal  Processing  Studies  and  Analysis- A  Study  of  electromagnetic  Inverse  Scattering. 
Syracuse  University  Research  Corporation,  Technical  Report  No.  RADC-TR-66-830,  April  1967,  AD  813  981. 

2lbkt. 
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where 


H$  is  the  scattered  magnetic  field  from  the  object. 

n  is  a  unit  vector  normal  to  the  surface  of  the  object. 

H  is  the  total  magnetic  field  at  the  surface  of  the  object. 

<t>  is  the  free  space  scalar  Green’s  function. 

Making  the  infinite  conductivity  assumption  for  the  object,  Bojarski  derives  the  follow¬ 
ing  result  for  r(z):3 


r(z)  = 


where 


—■  arg  P(6) 

“  °  '  -% 


cos'/2  0  d  9 ; 


(4) 


j  =v^" 

kQ  is  the  free  space  propagation  constant  equal  to  2ir/X,  where 
X  is  the  wavelength 

p  (0)  is  related  to  the  radar  cross  section  as  follows: 


P  P" 


(5) 


Norbert  N.  Bojarski.  Signal  Processing  Studies  and  Analysis-A  Study  of  Electromagnetic  Inverse  Scattering, 
Syracuse  University  Research  Corporation,  Technical  Report  No.  RADC-T R-66-830.  April  1967,  AD  813  981. 
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Since  p  is  complex,  the  amplitude  and  phase  of  the  received  signal  must  be  measured. 
The  parameter  p  is  called  the  complex  scalar  field  cross  section  and  can  be  related 
to  the  scattered  electric  field  as  follows: 


P 


2i/T  R 


(6) 


where 


R  is  the  far  field  distance  from  the  object  to  the  radar. 
E$  is  the  scattered  electric  field. 

Ej  is  the  electric  field  incident  on  the  object. 


The  three  equations  (1),  (2),  and  (4)  provide  two  methods  for  deten  'g  the 
profile  function  r(z)  of  an  infinitely  conducting  object  of  revolution. 


These  kinds  of  inverse  scattering  problems  appear  to  not  apply  to  the  problem  of 
terrain  feature  extraction,  because  terrain  consists  of  many  scattering  objects  of 
complicated  geometry  and  finite  conductivity. 


Profile  Inversion.  Many  inverse  problems  in  remote  sensing  can  be  formulated 
in  terms  of  a  Fredholm  integral  equation  of  the  first  kind  as  given  below: 


g(y) 


J  K(y,  x)  f  (x)  dx 
a 


(7) 
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In  the  above  equation,  K(y,  x)  forms  the  kernel  of  the  integral.  Measurements  are 
made  at  various  values  of  y  to  give  g.  The  objective  of  these  problems  is  to  determine 
f(x)  from  a  knowledge  of  K(y,  x)  and  a  finite  set  of  measurements  of  g(y).  Twomey 
has  written  about  techniques  for  the  solution  of  the  above  integral  equation.4  Even 
if  we  assume  perfect  precision  in  making  the  measurements  of  g(y),  the  questions 
of  existence  and  uniqueness  still  exist.  Equation  (7)  represents  a  mathematical  model 
of  some  physical  process  and,  as  such,  contains  simplifications  and  approximations  that 
may  or  may  not  be  valid  at  the  time  the  measurements  for  g(y)  are  made.  If  the 
approximations  are  not  valid,  the  existence  of  any  useful  solution  is  extremely  doubt¬ 
ful.  even  if  the  measurements  contain  no  errors.  If  the  model  is  shown  to  be  valid, 
then  one  must  address  the  question  of  uniqueness.  Given  a  set  of  measurements  for 
g(y).  a  solution  for  f(x>  can  be  determined,  but  the  question  is  whether  or  not  the 
solution  is  unique.  Other  solutions  may  exist  that  are  just  as  valid  and  physically 
realizable.  Many  times,  experimental  measurements  contain  errors  of  various  kinds, 
and  for  certain  problems  employing  equation  7,  this  leads  to  ill-posed  problems. 
An  ill-posed  problem  is  defined  as  one  in  which  small  measurement  errors  in  g(y) 
lead  to  large  fluctuations  in  f(x).  Although  the  difficulties  associated  with  existence, 
uniqueness,  and  ill-posed  problems  have  been  investigated  many  times,  they  are 
still  the  subjects  of  much  research. 


Example  of  a  Real  Profile  Inversion  Problem.  In  determining  earth  conduc¬ 
tivity  with  depth  surface  current  measurements  are  made  using  two  parallel  conductors 
lying  on  the  surface  of  the  ground  with  a  separation  distance  x  and  a  potential  dif¬ 
ference  v  applied  between  them.  The  resulting  electric  field  distribution  1  (x.z)  is  a 
function  of  the  separation  distance  x  and  the  depth  z  in  the  soil.  This  electric  field 
distribution  can  be  computed  if  the  permittivity  of  the  soil  is  constant  oi  varies  in  a 
known  way.  The  conductivity  of  the  soil  a(z)  van  s  with  depth  arid  can  be  related  to 
the  current  i(x)  in  the  conductors  as  follows: 


i(x) 


CO 


(8) 


S.  1  womey.  Introduction  to  The  Mathematies  of  Inversion  rn  K emote  Sensing  anil  Indireet  Measurements,  l  Iscvk’r 
Scientific  Publishing  Company.  1977. 
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By  making  measurements  of  the  current  and  the  electric  field  distribution  for  various 
values  of  x,  one  can  determine  the  conductivity  profile  if  the  integral  equation  can 
be  inverted. 


fo  obtain  a  numerical  solution  to  an  integral  equation  like  the  Fredholm  equa¬ 
tion,  one  must  rely  on  the  approximation  of  the  integral  by  a  sum,  that  is,  an  appro¬ 
priate  numerical  integration  formula,  such  as  Simpson’s  Rule,  must  be  chosen  and  the 
integral  equation  given  in  (7)  reduced  to  a  matrix  equation  of  the  following  form: 


g  =  A  x  f 


(9) 


where  g  is  an  Mxl  matrix  corresponding  to  the  M  measurements  that  have  been 
made.  The  objective  is  now  to  solve  the  matrix  equation  in  (9)  for  the  appropriate 
values  of  f. 


Calculating  of  the  Dielectric  Constant  for  a  Slightly  Rough  Surface.  Now 
consider  an  elementary  example  of  an  inverse  scattering  problem  in  which  a  solution 
for  the  dielectric  constant  of  a  slightly  rough  surface  is  obtained  in  terms  of  the  hori¬ 
zontally  and  vertically  polarized  radar  backscatter  coefficients.  The  scattering  geo¬ 
metry  of  the  problem  to  be  solved  is  shown  in  figure  2. 


A  plane  wave  with  a  harmonic  time  dependence  of  exp(ju>0  is  incident  onto 
the  slightly  rough  dielectric  surface  s(x.  y)  at  an  angle  of  incidence  (h.  The  medium 
above  the  surface  s(x.  y)  is  assumed  to  be  free  space  and  has  the  propagation  con¬ 
stant  ko.  The  propagation  constant  in  the  medium  below  the  surface  is  k.  The 
equations  for  the  radar  backscatter  coefficient  for  horizontal  and  vertical  polarizations 
are  given  as 

a°h(0.)  =  8o2k4  eos40.  R2  W(2kQ  sinfl.,0)  (10) 

<v<V  =  (ID 

“°l  k^cos^tk2  -  k2  )2  ( kcos0<  1  -  R„)  cosO,  +  k  sii^Od  I  +  RM)]2W(2k  sin0  0) 

k2  (k  cos  0.  +  k(i  cos  0)2 
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FIGURE  2.  Slightly  Rough  Surface. 

where 


o°h(0  )  is  the  radar  backscatter  coefficient  for  the  case  of  hori¬ 
zontal'  polarization  transmit  and  horizontal  polarization  receive, 
o°  (0.)  is  the  radar  backscatter  coefficient  for  the  case  of  vertical 
polarization  transmit  and  vertical  polarization  receive, 
of  is  the  standard  deviation  of  the  rough  surface  undulations  of 
s(x,  y) 

k  cos  0  =  >/k2  —  k2  sin2  0. 


W(2kc  sin0j,O)  is  the  surface  roughness  spectrum  and  is  equal 
to  the  two-dimensional  Fourier  transform  of  the  autocorrelation 
function  of  s(x,  y),  and  R±  and  R,,  are  the  Fresnel  reflection 
coefficients  for  horizontal  and  vertical  polarizations,  respectively, 
and  are  given  as 
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(12) 


k  cos  0 .  -  %/k2  -  k2  sin2  0 
k  cos  0.  +  x/k2  -  k*  sin2  6 . 

O  I  v  O  I 

k2  cosO.  —  k  y/v}  -  k2  sin2  6. 

I  O  v  O  I 

k2  cosO.  +  k  v/k2  -  k2  sin2  0. 

I  O  v  O  1 


(13) 


k  =  <jj  y/Ji  e  e~  where  e  is  the  relative  dielectric  constant  of 

v  *o  or  r 

the  medium  below  the  rough  surface  s(x,  y). 


Equations  (10)  and  (11)  are  derived  from  the  first-order  perturbation  theory, 
which  uses  a  two-dimensional  Fourier  transform  representation  for  the  scattered 
Fields.  The  basic  requirement  that  equations  (10)  and  (11)  must  obey  is  that  k0o, 
must  be  less  than  one.  This  is  because  a  perturbation  approach  is  used  for  the  solution, 
and  only  firs'-order  terms  are  considered.  In  equations  (10)  and  (11),  the  back- 
scatter  coefficients  are  not  only  functions  of  the  dielectric  constant  (er)  and  incidence 
angle  (0j)  but  also  functions  of  the  surface  roughness  parameters.  Any  attempt  to 
determine  the  dielectric  constant  from  either  (10)  or  (11)  will  require  a  knowledge 
of  the  backscatter  coefficient  and  of  the  surface  roughness  properties.  A  solution 
should  be  obtained  for  the  dielectric  constant  (ef)  in  terms  of  remote  sensing  mea¬ 
surements  such  as  the  backscatter  coefficient.  The  dielectric  constant  can  be  directly 
related  to  moisture  content,  and  if  a  solution  for  the  dielectric  constant  can  be  found 
in  terms  of  backscatter  coefficient,  a  quantitative  way  of  determining  soil  moisture 
will  have  been  found. 


However,  it  cannot  be  assumed  that  knowledge  of  the  rough  surface  properties 
exists  for  a  typical  terrain  surface  for  which  equations  (10)  and  (11)  are  valid. 
On  the  other  hand,  if  equation  (10)  is  divided  by  equation  (11),  the  rough  surface 
properties  cancel  out  and  are  no  longer  required. 


Let 


f.  (*•)  =  (0.)/a°  (0.) 

hv  i  hh  if  vv'  i 


(14) 
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then 


4k2  cos2  9  R2  k2  (k  cos  6  +  k  cos4>)2 

O  1  X  IQ  -  — 

(k2  -  k2)2[k  cos^  (1  -  R„)  cosOj  +  kQ  sin2  9.  (1  +  Rn)l 2 


In  figure  3,  the  results  are  presented  for  computing  fhv(0j)  for  two, values  of 
e  over  an  angular  range  of  0  to  80  degrees.  We  see  that  when  9t  =  0  ,  fhv(0j) 
is'  equal  to  one  for  both  dielectrics.  This  is  because  when  9{  is  equal  to  0°,  the 
horizontal  and  vertical  polarizations  will  give  the  same  result.  As  the  angle  of  incidence 
is  increased,  the  larger  value  of  er  makes  fhv(0j)  fall  off  faster.  Taking  the  square 
root  of  (15)  leads  to  the  following  result: 

f(k)  =  (k2  -k2)>/  fhv(flj)  [kcos<*>(l  -  R^cosOj  +  k2sin204(l  +  Ru)l 

-2k  R  cos0.  (k2cos0.  +  k  kcos<»  =  0  (16) 

Oil  I  O 


Equation  (16)  can  be  used  to  obtain  a  solution  for  k  and.  hence.  et  by 
employing  the  Newton-Raphson  numerical  methods  technique.  In  this  method,  the 
n  +  1st  estimate  of  k  is  related  to  the  n‘h  estimate  as  follows: 


k 


n  +  l 


f'<kn) 


(17) 


The  prime  on  f  is  used  to  indicate  a  derivative  of  f  with  respect  to  k.  An  expression 
for  this  derivative  can  be  written  as  follows: 


f(k)  =  _£!_  =  -2kyfhv(d7  (kcos0(l -R||)cos0.  + ko sin2  0.(1 +Rn)l 


dk 


+  /w<k!-k2)  p(k^{1_R  ^ 
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FIGURE  3.  Study  of  fhv(0j)  Variations  With 


1  r  dRx  2 

-  kcos^cosfl.)  -  2k  cost).  - Mk  cost).  +  k  k  cos0) 

1  J  0  •  L  dk  10 


d(k  cos4>) 

+  R,(2kcos0.  +  k  - ) 

1  1  0  dk 


'] 


dR  dR.. 

To  compute  T  (k),  there  is  a  need  to  calculate  — — — —  and 


(18) 

d(kcos$) 


dk  ’  dk  dk 

These  three  derivatives  can  be  computed  from  the  equations  given  previously. 


where 


dRt  _  -  2k  kQ  cos  0{ 

dk  (k2  -  k2sin20.),/j  jkocos0j +~^7k^-^k2sin20.  j 2 

'  09) 


dR„  _  2kc  k  cosdj  [2W2  -  k2  1 
dk  W  {k2  cos0s  +  ko  W)2 


(20) 


W  =  -J  k2  -  k2  sin2  0. 


d(k  cos 0) 
dk 


K  ~ 

-  -  -  -  —  —  ■  =  k/W 

y/  k2  -  k2  sin2  0. 

V  O  I 


(21) 


To  validate  using  (17)  to  compute  k  and  er,  a  computer  program  was  writ¬ 
ten  to  solve  the  direct  scattering  problem  from  equations  (10)  and  (11).  Parti¬ 
cular  values  were  chosen  for  <7,.  and  er  A  form  for  the  roughness  spectrum 

W(2kosin0j.0)  was  assumed  and  values  were  calculated  for  o°h  and  o“v.  These 
values  were  then  used  to  compute  fh v (0 from  (14).  This  value  of  fhv(0j)  was 
used  to  calculate  k  from  ( 1 7).  It  was  found  that  for  values  of  greater  than 
zero,  er  could  be  determined  from  (17).  Thus,  the  inversion  technique  and  equations 
were  verified.  If  further  work  on  this  problem  is  warranted,  experimental  measure¬ 
ments  should  be  conducted,  perhaps  in  a  desert  area  to  determine  o£h  and  o°y 
as  a  function  of  incidence  angle.  Also,  ground  truth  data  would  have  to  be  collected, 
such  as  dielectric  constant  and  soil  moisture  content.  This  approach  would  enable 
the  inversion  techniques  to  be  analyzed  using  experimental  data. 
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Calculating  the  Dielectric  Constant  and  Conductivity  for  a  Medium  Beneath 
a  Forest  or  Jungle  Canopy  —  Horizontal  Polarization.  Let’s  now  consider  an 
inverse  problem  in  which  a  solution  is  sought  for  the  dielectric  constant  and  the 
conductivity  of  a  medium  beneath  a  forest  or  jungle  canopy.  The  geometry  for  this 
problem  is  presented  in  figure  4. 


The  simplest  electromagnetic  model  for  a  forest  canopy  consists  of  a  lossy 
dielectric  slab  over  a  flat  earth.  The  forest  is  of  height  d  and  is  considered  to  have 
an  average  dielectric  constant  (e2)  and  an  average  conductivity  (o2).  This  simple 
model  for  a  forest  or  jungle  is  a  good  approximation  for  frequencies  less  than  200 
MHz.  Although  even  in  this  frequency  range,  one  must  often  incorporate  anisotropic 
effects.  The  propagation  of  electromagnetic  waves  in  forested  media  has  been  studied 
extensively  for  frequencies  less  than  200  MHz.  (see  Wait,  et  al.)s  These  studies 
were  conducted  to  further  the  effort  of  point-to-point  communications  in  a  forest 
environment.  Our  problem,  however,  is  to  invert  the  remote  sensing  equations  to 
obtain  a  solution  for  the  electrical  properties  of  the  ground  (e( ,  a, )  in  terms  of  the 
far  field  reflectivity.  One  of  the  difficulties  in  dealing  with  problems  of  forested  media 
is  how  to  determine  the  average  electrical  parameters  e2 ,  and  o2 .  Hagn  discusses 
various  techniques  that  can  be  used  to  determine  the  electrical  properties  of  forests  and 
presents  some  results  for  two  different  jungle  sites  in  Thailand.6  In  figure  4,  a  bistatic 
radar  setup  is  indicated  in  medium  3  to  measure  the  reflected  electric  field.  Below  the 
forest,  the  ground  surface  is  assumed  to  be  fiat  with  a  dielectric  constant  e,  and  a 
conductivity  <7, .  A  horizontally  polarized  plane  wave  with  a  harmonic  time  de¬ 
pendence  of  exp(jwt)  is  incident  onto  the  forest  medium  at  an  angle  of  incidence  tL. 
The  total  electric  field  in  medium  three  (t3  (  r  ))  can  be  written  as  follows: 

Lj  (  r  )  =  ayE0  ( e'jk°z  £OS  0  ■  +  Rx  ejk°*  cos  ° j  l  e'jk°x  sinei 

(  )  (22) 


J.  R.  Wait.  R.  II.  Oil.  and  T.  lelfer.  (I.ditors),  “Work  Shop  on  Radio  System  in  I -orcstcd  and  or  Vegetated 
I  nvironments,”  Technical  Report  No.  ACC--ACO-I-74.  institute  lor  Telecommunication  Sciences.  0.  S. 
Department  of  C  ommerce,  liouldcr.  Colorado. 

Ilagn,  Ci.ll.,  Electrical  Properties  of  I'oreslrd  Media  in  the  "Work  Shop  on  Radio  System  in  I  orcstcd  and  or 
Vegetated  I- nvironments.”  Technical  Report  No.  ACC-ACO-1--74.  Institute  for  Telecommunication  Sciences, 
l'-S.  Department  of  Commerce.  Moulder.  Colorado. 


17 


where 


ay  is  a  unit  vector  in  the  y  direction. 

hu  is  the  amplitude  of  the  incident  electric  field. 

Rj  is  the  reflection  coefficient. 


The  first  term  on  the  right  side  of  (22)  represents  the  incident  field.  The  second 
term  represents  the  reflected  field.  The  reflection  coefficient  can  be  written  in  terms 
of  the  impedances  of  the  three  media  involved. 


where 


R, 


z.,2)  -  z 

inh 


3  It 


Z<2)  +  Z 

inh  3h 


(23) 


Z' 2> 

inh 


Z  (Z  +Z  )  +  Z  (Z  -7  X  e‘i2k  2d  tosfl  2 

i  h  ^  l  h  +  2  h  ^2  It  I  h  ^2h)e _ 

z.  +  z  -  (Z  ~Z  )e'i2k2d  cos  9  2 

lh  Z'2h  l/L|h  ^h'6 


(24) 


'3  h 


k„cos  d. 


'2  It 


k2cos  02 


'I  h 


k ,  cos  0 , 


kQ  sin  =  k,  sin  0  j 


kx  cos  02 


The  parameters  kQ,  k2 ,  and  kt  represent  the  propagation  constants  in  air. 
in  the  forest  layer,  and  in  the  ground,  respectively.  Although  k2  and  k,  are  both 
complex,  kQ  will  be  real.  Because  we  arc  dealing  with  complex  propagation  constants, 
the  reflection  coefficient  will  also  be  complex  and  can  be  written  as 


R 


i 


(25) 
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where  pL  represents  the  magnitude  of  Ri  and  <j>L  is  the  phase.  In  Figures  5  through 
8,  curves  are  presented  of  pL  versus  incidence  angle  for  three  different  materials 
below  the  forest  canopy.  The  three  materials  are  loam  soil,  water,  and  aluminum. 
The  electrical  properties  of  these  materials  are  given  below: 


erl  =  el  /  e0 

o,  (mhos  /  meter) 

Loam  Soil 

13.77'/  Water 

14.5 

0.0105 

Water 

80 

0.000234 

Aluminum 

1.0 

3.53  X  107 

The  values  used  for  e1  and  a2  are  1.01e0  and  4  x  10  s  mho/meter,  res¬ 
pectively.  These  numbers  were  taken  from  Hagn  and  represent  the  parameters  for  a 
jungle  site  in  Thailand.7  The  height  of  the  forest  was  assumed  to  be  20  meters. 
Figure  5  presents  curves  for  a  frequency  of  30  MHz.  The  oscillations  that  occur 
in  the  curves  are  due  to  the  influence  of  the  forest  layer.  At  large  angles  of  incidence, 
the  three  curves  tend  to  merge,  indicating  that  the  three  materials  would  be  indis¬ 
tinguishable  from  one  another.  As  the  frequency  decreases,  notice  that  the  forest 
layer  has  less  and  less  influence  on  pi ,  until  finally  in  Figure  8,  at  3  MHz,  all  three 
curves  are  completely  monotonic.  Also,  in  Figure  8,  the  water  and  the  loam  soil  would 
be  indistinguishable  for  all  angles  of  incidence.  In  figure  9  and  10,  curves  are  shown 
of  0L  in  degrees  versus  angle  of  incidence  for  the  loam  soil  and  water.  In  figure  1 1 , 
curves  are  shown  of  pL  versus  frequency  for  the  three  materials  at  an  angle  of  in¬ 
cidence  of  30°.  In  attempting  to  obtain  a  solution  for  e,  and  ax,  assume  that  both 
the  amplitude  and  phase  of  the  reflected  field  are  measured.  This  would  then  enable 
a  complete  determination  to  be  made  of  Ri.  Using  this  knowledge  of  R^,  we  will 
seek  to  invert  equations  (23)  and  (24)  and  obtain  a  solution  for  e,  and  o,  in 
terms  of  Rr  d,  e2,  o,,  0t  and  frequency.  By  combining  equations  (23)  and 
(24),  one  obtains: 


llagn.  C.  H..  Electrical  Properties  of  forested  Media  in  the  "Work  Shop  on  Radio  System  in  f  orested  and  or 
Vegetated  I  nvironments,”  Technical  Report  No.  ACC-ACO-1-74.  Institute  for  Telecommunication  Sciences. 
U.S.  Department  of  Commerce.  Boulder.  Colorado. 
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Study  of  Reflection  Coefficient  Versus  Angle  of  Incidence  -  Horizontal  Polarization. 


FIGURE  6.  Study  of  Reflection  Coefficient  Versus  Angle  of  Incidence  -  Horizontal  Polarization 


Polarization  Horizontal 


FIGURE  8.  Study  of  Reflection  Coefficient  Versus  Angle  of  Incidence  --  Horizontal  Polarization. 


Polarization:  Horizontal 
f  *  30  MHz  d  -  20  meters 


FIGURE  9.  Study  of  Phase  Versus  Angle  of  Incidence  for  Loam  Soil  -  Horizontal  Polarization. 


FIGURE  11.  Study  of  Reflection  Coefficient  Versus  Frequency  -  Horizontal  Polarization. 


,(2)  -  Z3h  *  *  +  Rj  *  _  Z2h  ^Zl  h  +  Z2h^  *  Z2h  ^Z1  h  ~Z2h^  exP^~2jk2d  COS  fl2) 


in  h 


1  -  R, 


zih  +  z2h  "  <zih -z2h)exp(-2jkjd  cos#,) 


(26) 


Using  (26)  to  solve  lor  Zlh  produces 


Z 


1  h 


z'h " -R.> 


[- 


-2jk2<lcos  »2 


]- 


Z2h  Z3h 


+  Rx)l  1  +e 


WM0 


- 2 j k 2dcos  a  2 


°l  =  OJM0 


ZjhO  +  R,) 

j  e-2jk  2dcos  e  2 

-z2h(i-Rl) 

['♦ 

e~2 jk  2 d cos  0  2 

Z22h(l-RL) 

"l  _e-2Jk2dcos  ®2  " 

-Zjh^hf +Ri> 

'l  +  e-2jk 2d c°s  e  2j 

(30) 


where 


Re  means  take  the  real  part  of  the  expression  and 
lm  means  take  the  imaginary  part  of  the  expression. 


The  above  two  equations  represent  the  final  result  of  the  derivation.  The  key  to 
using  (29)  and  (30)  is  in  measuring  R  .  Experimental  work  would  have  to  focus 
on  how  well  the  amplitude  and  phase  of  the  reflected  field  could  be  measured  in 
determining  e,  and  a{  .  Equations  (29)  and  (30)  were  validated  by  choosing 
particular  values  for  frequency.  Oj,  o2,  d.  e,.  and  or  which  were  then  used 
to  calculate  Rx .  This  calculated  value  of  Rx  was  treated  as  having  been  obtained 
from  measurements  and  was  used  in  (29)  and  (30).  The  resulting  values  for  e, 
and  o,  were  then  compared  to  the  original  values  used  to  compute  R^  It  was  found 
that  equations  (29)  and  (30)  provide  a  valid  means  of  solving  for  e,  and  a{. 


Calculating  the  Dielectric  Constant  and  Conductivity  for  a  Medium  Beneath 

a  Forest  or  Jungle  Canopy  -  Vertical  Polarization.  Let's  consider  the  forest 

canopy  problem  for  the  case  where  the  incident  wave  is  vertically  polarized.  The  geo¬ 
metry  of  the  problem  remains  as  shown  in  figure  4.  The  total  magnetic  field  in 
medium  three  (  (  r  ))  can  be  written  as  follows: 


H3  (  r  ) 


=  a, 


"o  {' 


-  j  k  /.cos  II ; 


+  R, 


jk 


Ozcos  o 


-jk0ssin  ii  | 


(31! 
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where 


H0  is  the  amplitude  of  the  incident  magnetic  Held 
R|(  is  the  reflection  coefficient 


The  reflection  coefficient  R|(  can  be  written  in  terms  of  the  impedances  of 
the  three  media  involved : 


R„  “ 


Z,  -  Z(2) 

3  v  in  v 

z  +  z<2) 

3  v  in  v 


(32) 


where 


,<2>  = 


Z2v(Zlv+Z2y)  +  Z2v  (Z,v  -Z2v)  exp  (-2jk2dcos02) 
Ziv  +Z2v  “  (Ziv  -Z2v)exp(-2jk2dcos02) 


Zlv  = 


kucos  0{ 


'2  v 


k2cos  9 2 
w(e2  -jo2/w) 


-lv 


k,cos0j 
w(e,  -jo,/cj) 


kosin0j  =  kjsinfl,  =  k2sin02 


(33) 


Once  again,  propagation  constants  for  media  1  and  2  are  complex.  This  will 
make  the  reflection  coefficient  complex  so  that  it  can  be  written  as 


R, 


(34) 


where  pM  is  the  magnitude  of  R|(  and  <t>{]  is  the  phase.  In  figures  12  through 
I  7,  curves  of  pn  versus  angle  of  incidence  are  shown  for  the  three  materials  discussed 
previously.  The  influence  ot  the  lorest  layer  is  not  nearly  so  obvious  in  this  case  as 
it  was  lor  the  horizontal  polarization  case  where  the  oscillations  in  the  curves  clearly 
indicated  the  presence  of  a  layer.  The  dip  that  occurs  in  the  curves  between  the  angles 
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Polarization  Vertical 
f  »  JO  MHz  d  =  20  meters 


FIGURE 


Polarization:  Vertical 


FIGURE  13.  Study  of  Reflection  Coefficient  Versus  Angle  oflncidence  -  Vertical  Polarization 


FIGURE  14.  Study  of  Reflection  Coefficient  Versus  Angle  of  Incidence  -  Vertical  Polarization. 


FIGURE  1 5.  Study  of  Reflection  Coefficient  Versus  Angle  of  Incidence  Vertical  Polarization. 


FIGURE  16.  Study  of  Reflection  Coefficient  Versus  Angle  of  Incidence  Vertical  Polarization. 


FIGURE  1 7.  Study  of  Reflection  Coefficient  Versus  Angle  of  Incidence  -  Vertical  Polarization. 


of  70°  and  80°  is  a  Brewster  angle  effect,  in  figures  18  and  19,  carves  of  phase 
versus  angle  of  incidence  are  shown  for  loam  soil  and  water.  In  figure  20,  curves  of 
p  versus  frequency  are  shown  for  the  three  different  materials  and  at  an  angle  of 
incidence  of  30  degrees.  As  in  the  horizontal  polarization  cases,  assume  that  the 
magnitude  and  the  phase  of  the  reflection  cofficient  are  measured  and  are  available 
to  determine  the  electrical  properties  e,  and  or  Using  equations  (32)  and  (33), 
one  obtains  an  expression  for  Z ,  v : 


'1  V 


Z,V[Z3V  (1  -Rn)  (1  +  e~2ik2dcos  0  2 )  -  z2v  (1  +  R„)  (1  -c'2jk2dcos92)] 
z2v  (1  +  Rh)  (1  +  e'2ik2dcos  0  2  >  _z3v(l  -Rn)  (1  _e'2jk2dcos  0  2 ) 


(35) 


Letting  the  quantity  on  the  right  side  of  (35)  be  equal  to  g  and  using  the  definition 
of  Z.  results  in 


k,  cos  0 , 
oat,  -  jo, 


g 


(36) 


Multiplying  the  numerator  and  denominator  on  the  left  hand  side  of  (36)  by  top0 
and  realizing  that  k2  =  +  oa2M  )eI  one  obtains 


cap  cos  0 . 

g  =  - - - -  (37) 

k. 

Squaring  (37).  multiplying  the  numerator  and  denominator  on  the  right  side  of  (37) 
by  k2 ,  and  substituting  1  -  sin2  0[  for  cos2  0,  results  in 

ca2p2  k2  { 1  -  sin2  0  ) 

_  O  1 _ I 


(38) 
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Polarization:  Vertical 


FIGURE  19.  Study  of  Phase  \  ersus  Angle  of  Incidence  for  Water  -  Vertical  Polarization 


FIGURE  20.  Study  of  Reflection  Coefficient  Versus  Frequency  -  Vertical  Polarization. 


r 


Recognizing  that  k,  sin  0,  =  ku  sin  0r  one  can  write  equation  (38)  as  follows: 


.42  22.2,  22.2.  2  a  „ 

k.g  -  o)  n  k.  +  co  n  k  sin  0  -  0 

1  D  O  l  'o  O  1 


(39) 


Equation  (39)  is  a  fourth  degree  equation  for  k. .  which  can  be  solved  by  letting 

V  * 


k;  = 


j-2  2  22v.  2  2  ,  2  .  2  ~  _  n 

£  g  -  w  p  {  +  w  |i  k  Sin  0.  =  0 

s  e  o  s  o  o  i 


(40) 


A  solution  for  £  can  now  be  obtained: 


£  = 


or  /r  ±  y/ t 0*nA0  -  4g2  (O2  n\  k2  sin2  0~ 

2g2 


(41) 


The  problem  of  which  sign  to  use  in  (41)  is  easily  resolved  by  seeing  what 
happens  at  0j  =  0°.  When  0i  =  0°  in  (41)  and  the  minus  sign  is  used,  the  totally 
unrealizable  result  of  £  =  0  is  obtained.  This  indicates  clearly  that  the  proper  sign  to 


choose  in  (41)  is  the  plus  sign.  Because  £  =  k2  we  now  have 


-jwp0ci,  -t  w2M0e,  =  w2m2 


1  +  v/l-4g2  eosin2  OJ  n0 


2g“ 


(42) 


In  the  above  equation.  kn  =  u >  v  M0  e<> •  *s  now  possible’  to  solve  for  e, 


and  a 


l  • 


el  =  fcrl  eo  = 


1  +  J  1  -4g2  eusin2  (T  j  fi~ 


28- 


(43) 
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1  +  >/ 1  —  4g2  eGsin2  0{  /  p0 


(44) 


u,  =  -top  1m 

1  *o 


Equations  (43)  and  (44)  are  die  final  desired  results.  One  can  see  how  for¬ 
midable  the  results  are  if  the  right  hand  side  of  (35)  is  substituted  for  g.  The  above 
equations  were  validated  using  the  same  method  that  was  used  for  horizontal  polari¬ 
zation.  Particular  values  were  chosen  for  frequency,  0 ^  e2,  ,  d,  e,,  and  o, ; 

Then,  these  were  used  to  compute  R(  .  The  computed  value  of  R(|  was  then  con¬ 
sidered  as  a  measured  quantity  and  used  in  (43)  and  (44)  to  compute  6,  and  o, . 
These  resulting  values  for  e,  and  Oj  were  compared  to  the  original  values  of  e, 
and  a,  that  were  used  to  compute  R(|.  By  following  this  method,  it  was  shown  that 
(43)  and  (44)  are  a  valid  way  of  computing  e,  and  o,. 


In  the  following  section,  the  results  that  were  obtained  will  be  discussed  along 
with  using  radar  imagery  to  extract  terrain  feature  information. 


DISCUSSION 


In  the  previous  section,  several  examples  of  inverse  problems  were  considered, 
f  irst,  inverse  scattering  was  examined  for  discrete  objects.  However,  this  particular 
problem  has  very  little  application  toward  determining  terrain  feature  parameters 
because  terrain  consists  of  collections  of  scatterers  that  are  placed  more  or  less  ran¬ 
domly  in  space.  The  physics  of  the  problem  of  scattering  from  the  collection  of  objects 
is  different  than  that  for  a  single  object  by  itself.  Also,  the  infinite  conductivity  as¬ 
sumption  is  not  valid  for  any  natural  terrain  feature. 


Next,  profile  inversion  was  examined.  This  type  of  inverse  problem  also  appears 
to  have  little  application  in  determining  terrain  feature  parameters,  unless  one  is 
interested  in  computing  the  dielectric  or  conductivity  profile  in  the  earth’s  surface. 
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In  addition,  the  problem  of  calculating  the  dielectric  constant  of  a  slightly  rough 
surface  using  radar  backscatter  data  from  two  polarizations  was  considered.  The 
solution  for  this  problem  might  be  useful  in  determining  the  moisture  content  in  a 
soil.  An  experimental  measurements  program  would  have  to  be  set  up  to  see  how  well 
the  dielectric  constant  and  moisture  content  could  be  determined  from  actual  data. 
Also,  the  measurements  program  would  have  to  include  the  collection  of  ground 
truth  data  in  the  way  of  dielectric  constants,  moisture  contents,  and  perhaps  surface 
roughness  properties.  The  initial  experimental  work  could  be  performed  using  a  radar 
system  such  as  the  one  at  the  University  of  Kansas. 


The  final  problem  considered  was  the  calculation  of  the  electrical  properties 
beneath  a  forest  or  jungle  canopy.  A  solution  for  this  problem  is  important  when 
one  is  trying  to  ascertain  the  nature  of  terrain  materials  located  under  a  forest  canopy. 
The  electrical  properties  of  the  ground  were  computed  in  terms  of  a  measured  re¬ 
flection  coefficient,  the  height  of  the  forest,  and  the  electrical  properties  of  the  forest. 
Once  again,  a  measurements  program  would  have  to  be  conducted  in  order  to  test 
the  technique.  One  type  of  inverse  problem  that  has  not  been  considered  as  yet  is  that 
of  using  side-looking  airborne  radar  imagery  to  infer  various  terrain  features.  This 
problem  will  be  considered  in  the  following  paragraphs. 


The  average  gray  shade  on  a  radar  image  over  an  area  the  size  of  a  resolution 
cell  is  related  to  the  value  of  the  radar  backscatter  coefficient  (  a°  ).  This  backscatter 
coefficient  is  a  function  of  both  system  and  terrain  parameters  as  indicated  below: 


a°  =  f(x,0„p.,pr,rs,vs,su,a,zx,zy) 


Theoretical  expressions  for  a°  have  been  derived  for  certain  types  of  rough 
surfaces,  vegetation,  and  snow.  The  theoretical  efforts,  however,  are  just  beginning  to 
explain  certain  aspects  of  measured  data,  and  much  effort  is  needed  in  this  area  yet. 
In  the  above  expression  for  a°,  the  following  notation  has  been  used. 
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Radar  System  Parameters 


X  =  radar  wavelength 

6.  =  local  angle  of  incidence 

Pj  =  polarization  of  the  incident  wave 

Pf  =  polarization  of  the  received  wave 


Terrain  Parameters 

r  =  surface  roughness  properties 

Vs  =  volume  scattering  affects  such  as  those 
occurring  in  vegetation  and  snow. 

Su  =  subsurface  scattering  affects 

e  =  complex  dielectric  constant  of  the  terrain  surface 
and  Zy  =  two  orthogonal  slopes  of  the  terrain  surface 

Variations  in  gray  shade  from  one  resolution  cell  to  another  are  due  to  changes 
in  o° ,  which  in  turn  is  due  to  changes  in  one  or  more  of  the  above  parameters.  Inter¬ 
pretating  radar  imagery  consists  of  trying  to  infer  terrain  features  from  the  tonal 
variations  on  the  image.  The  radar  image  interpreter  should  have  a  good  understanding 
of  how  each  of  the  system  and  terrain  parameters  influences  a° .  This  understanding 
can  be  obtained  by  analyzing  and  developing  theoretical  models  and  by  analyzing 
experimental  data.  A  great  deal  of  research  effort  has  gone  into  radar  image  inter¬ 
pretation  using  human  photointerpreters.  Also,  people  who  were  trained  as  interpreters 
of  aerial  photography  have  been  used  to  analyze  radar  imagery.  However,  these  photo¬ 
interpreters  do  not  have  the  necessary  theoretical  background  to  enable  them  to  pin¬ 
point  specific  parameter  influences,  which  has  tended  to  keep  radar  image  interpre¬ 
tation  largely  in  a  heuristic  state.  One  area  of  radar  image  analysis  that  has  not  received 
much  attention  is  that  of  applying  digital  image  processing  methods.  Analyzing  radar 
imagery  with  various  texture  measures  appears  to  be  worthy  research,  if  appropriate 
ground  truth  information  can  be  obtained. 
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CONCLUSIONS 


As  a  result  of  this  study,  the  following  conclusions  can  be  made: 

1.  Inverse  scattering  solutions  for  discrete  objects  are  not  applicable 
in  determining  terrain  features. 

2.  Profile  inversion  solutions  appear  to  have  only  a  slight  applicability 
in  determining  terrain  features. 

3.  The  problem  of  determining  the  dielectric  constant  of  a  slightly  rough 
surface  can  be  solved  using  the  horizontal  and  vertically  polarized 
backscatter  coefficients,  along  with  a  Newton-Raphson  numerical 
methods  technique. 

4.  A  solution  for  the  dielectric  constant  and  the  conductivity  of  a 
medium  located  beneath  a  forest  canopy  can  be  obtained  in  terms 
of  the  amplitude  and  phase  of  the  reflected  wave. 
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